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Abstract 

Events with a final state consisting of two or more photons and large missing 
transverse energy have been observed in e + e~ collisions at centre-of-mass energies in 
the range 192 — 209 GeV using the OPAL detector at LEP. Cross-section measure- 
ments are performed within the kinematic acceptance of the selection and compared 
with the expectations from the Standard Model process e + e~ — > 1/1/77(7). No evi- 
dence for new physics contributions to this final state is observed. Upper limits on 
cr(e+e- -» XX) • BR 2 (X -» Y7) are derived for the case of stable and invisible Y. 
In the case of massive Y the combined limits obtained from all the data range from 
10 fb to 60 fb, while for the special case of massless Y the range is 20 fb to 40 fb. The 
limits apply to pair production of excited neutrinos (X = u*,Y = u), to neutralino 
production (X = x[>,Y = Xi) an d to super symmetric models in which X = Xi an d 
Y = G is a light gravitino. 
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1 Introduction 



We describe measurements and searches performed using a data sample of photonic events 
with large missing transverse energy collected with the OPAL detector in 1999 and 2000, 
the final two years of LEP operation. The events result from e + e~ collisions in the centre- 
of-mass energy range of about 192 — 209 GeV with a combined integrated luminosity 
of 426.5 pb -1 . When deriving cross-section limits on new physics processes, these data 
are combined with previously published data [1] taken at 189 GeV and corresponding to 
177.3 pb -1 . The present paper builds on past publications based on data samples collected 
at lower centre-of-mass energies [1-3]. The new data samples, taken at the highest energies 
achieved by LEP, provide discovery potential in a new kinematic regime with a large 
increase in integrated luminosity. Similar searches have been made by the other LEP 
collaborations [4]. 

The analysis presented here is designed to select events with two photons and significant 
missing transverse energy in the final state, indicating the presence of at least one neutrino- 
like invisible particle which interacts only weakly with matter. The event selection for this 
search topology is identical to that used in our most recent publication [1]. Within the 
Standard Model, such events are expected from the e + e~ — * 2/1/77(7) process. The selection 
is designed to retain acceptance for events with an additional photon, provided that the 
system formed by the three photons is consistent with the presence of significant missing 
transverse energy. 

This final-state topology is also sensitive to several new physics scenarios. In the context 
of the search for new physics, the emphasis in this publication is on general searches 
applicable to a broad class of models. To this end, a generic classification is used: e + e" — > 
XX where X is neutral and can decay radiatively (X — > Y7) and Y is stable and only 
weakly interacting. The limits presented for this generic process are applicable to a variety 
of physics searches. For the general case of massive X and Y this includes conventional 
supersymmetric processes (X = x°) Y = Xi)- There is particularly good sensitivity for the 
special case of M Y ~ 0. This is applicable both to the production of excited neutrinos 
(X = u*,Y = is) and to supersymmetric models in which the lightest supersymmetric 
particle (LSP) is a light gravitino and Xi is the next-to-lightest supersymmetric particle 
(NLSP) which decays to a gravitino and a photon (X = x?,Y = G). In the latter case, 
we also set limits on an example light-gravitino model [5]. The neutralino lifetime in such 
models is a free parameter. In this paper we address only the case of promptly decaying X. 

This search topology also has sensitivity to the production of two particles, one invisible, 
or with an invisible decay mode, and the other decaying into two photons. Such events 
might arise from the production of a Higgs-like scalar particle, S° : e + e~ — > Z°S°, followed 
by S°^» 77, Z° — > vv. The results of an OPAL search for this process, including the 
hadronic and leptonic Z° decays, have been separately reported [6]. Finally, this search 
topology can also probe WW77 quartic couplings in the e + e~ — > f e Pe77 process. The 
OPAL quartic gauge coupling measurements are described in [7]. 

This paper first describes the OPAL detector and the Monte Carlo samples used. A brief 
summary of the event selection will then be given, followed by cross-section measurements 
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and comparisons with Standard Model expectations. The new physics search results will 
then be discussed. 



2 OPAL Detector and Monte Carlo Samples 

The OPAL detector, which is described in detail in [8], contained a silicon micro- vertex 
detector surrounded by a pressurized central tracking system operating inside a solenoid 
with a magnetic field of 0.435 T. The barrel and endcap regions of the detector were instru- 
mented with scintillation counters, presamplers and a lead-glass electromagnetic calorime- 
ter (ECAL). The magnet return yoke was instrumented for hadron calorimetry and was 
surrounded by muon chambers. Electromagnetic calorimeters close to the beam axis mea- 
sured luminosity and completed the acceptance. 

The measurements presented here are based mainly on the observation of clusters of 
energy deposited in the lead-glass electromagnetic calorimeter. This consisted of an ar- 
ray of 9,440 lead-glass blocks in the barrel region, | cos 6*| < 0.82, with a quasi-pointing 
geometry and two endcap arrays, each of 1,132 lead-glass blocks, covering the polar an- 
gle 1 range, 0.81 < |cos#| < 0.984. Hermetic electromagnetic calorimeter coverage was 
achieved beyond the end of the ECAL down to 33 mrad in polar angle with the use of the 
gamma-catcher calorimeter, the forward calorimeter and the silicon-tungsten calorimeter. 

Scintillators in the barrel and endcap regions were used to reject backgrounds from 
cosmic-ray interactions by providing time measurements for the large fraction 80%) of 
photons which converted in the material in front of the ECAL. The barrel time-of-flight 
(TOF) scintillator bars were located outside the solenoid in front of the barrel ECAL 
and matched its geometrical acceptance |cos#| < 0.82. Tile endcap (TE) scintillator 
arrays were located in front of the endcap ECAL at 0.81 < |cos#| < 0.955. Additional 
scintillating-tile arrays, referred to as the MIP plug, were located at more forward angles. 
In the region from 125 to 200 mrad these detectors were used to provide redundancy in 
the rejection of events with significant electromagnetic activity in the forward region. 

The integrated luminosities of the data samples are determined to better than 1% from 
small-angle Bhabha scattering events in the silicon-tungsten calorimeter. Triggers based on 
electromagnetic energy deposits in either the barrel or endcap electromagnetic calorimeters 
lead to full trigger efficiency for photonic events passing the event selection criteria used 
in this analysis. 

The NUNUGPV98 [9] and KK2f [10] Monte Carlo generators were used to simu- 
late the Standard Model signal process, e + e~ — > 14777(7). For other expected Stan- 
dard Model processes, a number of different generators were used: RADCOR [11] for 
e+e" -> 77(7); BHWIDE [12] and TEEGG [13] for e+e" -> e + e-( 7 ); KORALW [14] 
using grc4f [15] matrix elements for e + e~ — > uu£ + £~('~f) and e + e~ — > ^z7qq(7), and KO- 
RALZ [16] for e + e~ — > fi + fi~('~f) and e + e~ — > t + t~ (7). The BDK program [17] was used 

1 The OPAL right-handed coordinate system is defined such that the origin is at the centre of the 
detector and the z axis points along the direction of the e~ beam. The polar angle 6 is defined with 
respect to the e~ beam direction and <f> is the azimuthal angle measured from the +x axis. 



5 



for e + e~ — > e + e~£ + £~, except for e + e~ — > e + e~e + e~ which was generated using the Ver- 
maseren program [18]. The expected contribution from each of these Standard Model 
processes was evaluated using a total equivalent integrated luminosity at least five times 
larger than the integrated luminosity of the data sample. 

To simulate possible new physics processes of the type e + e~ — > XX where X decays 
to Y7 and Y escapes detection, a modified version of the SUSYGEN [19] Monte Carlo 
generator was used to produce neutralino pair events of the type e + e~ — > X2X2! X2 ~^ X?7> 
with isotropic angular distributions for the production and decay of X2 an d including the 
effects of initial-state radiation. For = 206 GeV, Monte Carlo events were generated 
at 49 points in the kinematically accessible region of the (Mx, My) plane. Monte Carlo 
events at 42 points in (M x , M Y ) with y/s = 189 GeV were generated for our previous 
publication [1]. Using these two samples, the selection efficiency was determined for each 
generated point and then parametrized as a function of (M x , M Y ) and centre-of-mass 
energy. The efficiency varies slowly with energy and for energies above 206 GeV, the 
206 GeV values were used. All Monte Carlo samples described above were processed 
through the full OPAL detector simulation [20]. 



3 Event Selection 

A detailed description of the event selection is given in our previous publications [1,2]. In 
brief, photons are identified as energy deposits in the electromagnetic calorimeter. Events 
are required to have no other significant activity, except for the possibility of additional 
photons. Information from the tracking chambers is used to reject electromagnetic clus- 
ters associated with prompt charged tracks while retaining sensitivity for photons which 
converted in the material between the interaction point and the calorimeter. Timing infor- 
mation is used to reject backgrounds from cosmic-ray events. Events with activity beyond 
the acceptance of the ECAL are vetoed using information from the gamma catcher, the 
forward calorimeter, the silicon-tungsten calorimeter and the MIP plug. The kinematic ac- 
ceptance of the selection is defined by requiring: 

• at least two photons, each with x 7 > 0.05 and 15° < < 165°, or one photon with 
E 1 > 1.75 GeV and |cos#| < 0.8 and a second photon with E 1 > 1.75 GeV and 
15° < 6 < 165°; here E 1 is the photon energy, 6 is the photon polar angle and x 7 is 
the photon scaled energy E 7 /E heam 

• that the two-photon system consisting of the two highest-energy photons have mo- 
mentum transverse to the beamline (p 7 - 7 ) satifying j E beam > 0.05 

The selection is designed to retain acceptance for events with additional photons in 
which the resulting photonic system is still consistent with the presence of significant 
missing energy. This reduces the sensitivity of the measurement to the modelling of higher- 
order contributions. 
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4 Selection Results 



The data described in this paper were taken during the final two years of LEP operation, at 
centre-of-mass energies between 192 and 209 GeV. For the purposes of this publication the 
data have been binned into six samples with mean centre-of-mass energies of approximately 
192, 196, 200, 202, 205 and 207 GeV. The energy ranges and luminosity breakdown are 
summarized in Table Applied to the entire sample, the selection yields a total of 54 
events, in good agreement with the KK2f prediction of 57.2 ±1.3 events for the Standard 
Model e + e~ — > 1/7777(7) contribution. The expected contribution from other Standard 
Model processes and from cosmic ray and beam-related backgrounds is 1.2 ± 0.3 events, 
dominated by contributions from low-angle radiative Bhabha events and radiative four- 
fermion final states. The selection results are included in Tabled The selection efficiency 
for e + e~ — > ^7777(7) events within the kinematic acceptance of the selection is (65.7±1.5)%, 
independent of energy. The cross-section within the kinematic acceptance of the selection 
is also shown in Table ^ as are the corresponding predictions obtained using the KK2f 
Monte Carlo generator. The predictions of the NUNUGPV98 Monte Carlo generator were 
also examined and agreed well with those of KK2f. Small differences are accounted for in 
the systematic uncertainties. 

The dominant sources of systematic uncertainties arise from modelling of the event 
selection efficiency, especially the simulation of the detector material and consequent pho- 
ton conversion probabilities. The effects of these uncertainties and of uncertainties on the 
efficiency of timing cuts used to suppress cosmic-ray events are calculated accounting for 
different event topologies (both photons in the barrel region, both in the endcap, or one 
in each). This total uncertainty is 1.7%. Other sources arise from uncertainties on the in- 
tegrated luminosity measurement (0.5%), on detector occupancy estimates (1%) obtained 
from the analysis of randomly triggered events, on comparisons of different Monte Carlo 
event generators for the process e + e~ — > 1/7777(7) (1%). The total systematic uncertainty 
common to each energy bin is 2.3%. In individual energy bins, Monte Carlo statistics 
account for an additional systematic uncertainty of 0.9 — 1.4%. 

The kinematic properties of the selected events, summed over all energies, are displayed 
in Figure where they are compared with the predicted distributions for e + e~ — > 2/7777(7) 
obtained using the KK2f generator normalized to the integrated luminosity of the data. 
Plot (a) shows the recoil mass distribution of the selected events (for the two most energetic 
photons in the case of events with three or more photons). The distribution is peaked near 
the mass of the Z° as is expected for contributions from e + e~ — > 1/7777(7). The resolution of 
the recoil mass is typically 4 — 6 GeV for M reco ii ~ Mz- Events with a negative recoil- mass 
squared are plotted in the zero bin of the distribution. Plot (b) shows the distribution of 
the scaled energy of the second most energetic photon. Plot (c) shows the 77 invariant- 
mass distribution for which the mass resolution is typically 1 — 2 GeV. Plot (d) shows the 
distribution in scaled transverse momentum of the selected two-photon system. 

There are 3 selected events having a third photon with deposited energy above 300 MeV 
and within the polar-angle acceptance of the selection. The corresponding expectation from 
KK2f is 3.36 ± 0.08 events. 
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5 Data Interpretation 



The results of this selection are used to test the Standard Model and to search for new 
physics contributions. In the absence of an excess of events beyond the Standard Model 
expectation, we set 95% CL upper limits on the quantity cx(e + e~ — » XX) ■ BR 2 (X — > Y7) 
for the general case of massive X and Y, and separately for the special case of My ~ 0. 
Efficiencies were evaluated under the assumption that X decays promptly. Monte Carlo 
samples were generated for a variety of mass points in the kinematically accessible region 
of the (Mx, My) plane. To set limits for arbitrary Mx and My, the efficiency over the 
entire (Mx, M Y ) plane was parameterized using the efficiencies calculated at the generated 
mass points. For Mx values below Mz/2, search results based on LEP1 data have been 
previously reported [21]. In this low-mass region, events with radiative return to the Z° 
followed by Z° — > XX would yield very different kinematics than those used here to generate 
the signal Monte Carlo samples. For this reason, the search is restricted to the mass region 
M x > M z /2. 



5.1 Search for e + e — ► XX, X — ► Y7 ; General case: My > 

The searches for e + e~ — > XX, X — > Y7, both for the general case discussed here and 
the special case of M Y « discussed in section 5.2, use the methods described in our 
previous publications [1,2]. Selected events are classified as consistent with a given value 
of Mx and My if the energy of each of the photons falls within the region kinematically 
accessible to photons from the process e + e~ — > XX, X — > Y7, including resolution effects. 
Selection efficiencies at some of the generated grid points for the e + e~ — > XX, X — > Y7 
t/s = 206 GeV Monte Carlo events are shown in TableEl These values include the efficiency 
of the kinematic consistency requirement which is higher than 95% at each generated point 
in the region of the (Mx, My) plane. For Mx — My values lower than 5 GeV the efficiency 
begins to fall off rapidly and is thus difficult to model accurately. For this reason, we place 
limits only in the region of the (Mx, My) plane satisfying Mx — My > 5 GeV. Efficiencies at 
lower centre-of-mass energies are obtained from an interpolation between these efficiencies 
and the equivalent efficiencies at 189 GeV, which are given in our previous publication [1] . 
For data taken at centre-of-mass energies above 206 GeV, the 206 GeV efficiencies are used. 

Events from e + e~ — > 1/1/77(7) are typically characterized by a high-energy photon 
from the radiative return to the Z° and a second lower energy photon. The kinematic 
consistency requirement is such that the two photons must have energies within the same 
(kinematically accessible) region. Thus, as Mx and My increase, the allowed range of 
energy for the photons narrows, and fewer 1/1/77(7) events will be accepted. For the 54 
selected events, the distribution of the number of events consistent with a given mass 
point (Mx,My) is consistent with the expectation from e + e~ — > ^77(7) Monte Carlo, 
over the full (M x ,My) plane. Upper limits are placed on <r(e + e~ — > XX) • BR 2 (X — > Y7) 
accounting for the number of selected events and the expected number of background 
events from the process e + e~ — > ^7777(7). Other backgrounds are not subtracted. For 
each of the energy bins, Table H3 shows the maximum and minimum limits obtained in the 
region of the (Mx,My) plane described above. Figure El shows the 95% CL lower limits 
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on a(e+e- -»■ XX) ■ BR 2 (X -> Y7) at ^ = 207 GeV, obtained from all OPAL data with 
\fs > 189 GeV, under the assumption that cr(e + e~ — > XX) scales with centre-of-mass 
energy as flx/s. These limits range from 10 — 60 fb. 

Systematic uncertainties arise from the sources described in section 4. However there 
are additional contributions due to limited Monte Carlo statistics at each of the gener- 
ated (Mx, My) points and from uncertainties on the efficiency parameterization across the 
(Mx, My) plane and as a function of energy. The combined relative uncertainty on the 
efficiency varies from about 3% to 6% across the plane (for Mx — My > 5 GeV). The un- 
certainty on the expected SM background contribution is 2.6%. In calculating the limits, 
systematic uncertainties are accounted for in the manner advocated in reference [22]. This 
also applies to the limits for the My « case, presented in the next section. 

5.2 Search for e + e~ — ► XX, X — ► Y7 ; Special case: M Y « 

For the special case of My « the applied kinematic consistency requirements differ from 
those used for the general case. One can calculate [23] the maximum mass, M^, which is 
consistent with the measured three-momenta of the two photons, assuming a massless Y. 
A cut on Mx ax provides further suppression of the 1/1/77(7) background while retaining 
high efficiency for the signal hypothesis. This is discussed in more detail in reference [3]. 
To allow for resolution effects, we require that the maximum kinematically allowed mass 
be greater than Mx — 5 GeV. This has better than 96% relative efficiency for signal at all 
values of Mx while suppressing much of the remaining 1^77(7) background. 

The M™ distributions for all selected events, divided into the 192 — 202 GeV and 
205 — 207 GeV data samples, are shown in Figured In each case, the points with error bars 
show the OPAL data while the unshaded histogram shows the expected contribution from 
the e + e~ — > 1/1/77(7), from KK2f Monte Carlo, normalized to the luminosity of the data. 
Shown as a shaded histgram in the 205 — 207 GeV plot is the expected distribution from 
signal Monte Carlo events generated with Mx = 100 GeV (with arbitrary normalization). 
For this My « case, the signal reconstruction efficiencies calculated from Monte Carlo 
events generated at y/s = 206 GeV are shown in Table |U after application of the event 
selection criteria and then after the cut on M™ x - Also shown in Table |U are the numbers 
of events selected from the 205 — 207 GeV data sample which are consistent with each 
value of Mx as well as the expected number of e + e~ — > 2/7/77(7) events. The number of 
selected events (from the 205 — 207 GeV sample) consistent with a given value of Mx varies 
from 10, for Mx > 45 GeV, to 2 at the kinematic limit. The expected number of events 
decreases from 14.9 ±0.4 at M x > 45 GeV to 1.28 ±0.08 consistent with M x > 102.5 GeV. 

Based on the efficiencies and the number of selected events, we calculate 95% CL upper 
limits on <r(e + e~ — > XX) • BR 2 (X — > Y7) for My « as a function of Mx, in each 
region of centre-of-mass energy. The last two columns of Table |3] show the range of limits 
obtained from each of the data samples, for Mx values from 45 GeV up to the kinematic 
limit. Figure 0] shows the limit obtained from the 207 GeV data sample, as well as the 
combined limit obtained from the entire data sample with yfs > 189 GeV assuming that 
the cross-section scales as /3x/ s. For the mass range of interest (Mx > 45 GeV) the model- 
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independent limits range between 45 and 70 fb while the combined limits range between 20 
and 45 GeV. These limits 2 can be used to set model-dependent limits on the mass of the 
lightest neutralino in supersymmetric models in which the NLSP is the lightest neutralino 
and the LSP is a light gravitino (X = x^Y = G). Shown in Figure EJ as a dotted line, 
is the (Born-level) cross-section prediction from a specific light gravitino LSP model [5] in 
which the neutralino composition is purely bino, with mg H = 1.35m^o and mg L = 2.7m^o. 
Within the framework of this model, Xi masses between 45 and 99.0 GeV are excluded at 
95% CL. 

As described in section 2, the efficiencies over the full angular range have been obtained 
using isotropic angular distributions for the production and decay of X. The validity of 
this model has been examined based on the angular distributions calculated for photino 
pair production in reference [24]. For models proposed in reference [25], the production 
angular distributions are more central and so this procedure is conservative. For a 1 + cos 2 9 
production angular distribution expected for t-channel exchange of a very heavy particle 
according to reference [24], the relative efficiency reduction would be less than 2% at all 
points in the (Mx, My) plane. 



6 Conclusions 

We have searched for events with a final state consisting of two or three photons and large 
missing energy, in data taken with the OPAL detector at LEP, at centre-of-mass energies 
in the range of 192 - 209 GeV. 

The selection requires at least two photons with scaled energy x 1 > 0.05 within the 
polar angle region 15° < 9 < 165° or at least two photons with energy E 1 > 1.75 GeV 
with one satisfying Icos^l < 0.8 and the other satisfying 15° < 9 < 165°. In each case, 
the requirement p?P /E beiim > 0.05 is also applied. There are 54 events selected. The KK2f 
prediction for the contribution from e + e~ — > 1/1/77(7) is 57.2 ± 1.3 events; expected con- 
tributions from other sources sum to 1.2 ± 0.3 events. The number of events observed 
in the data and their kinematic distributions are consistent with Standard Model expec- 
tations. Limits on new physics processes of the form cr(e + e~ — > XX) ■ BR 2 (X — * Y7) are 
set separately at energies of 192, 196, 200, 202, 205 and 207 GeV. In addition, combined 
limits are set at \fs = 207 GeV, assuming a /3%/s scaling of the production cross-section 
cx(e+e- -> XX). From the full OPAL data sample with y/s > 189 GeV, we derive 95% CL 
upper limits on <r(e + e~ — > XX) • BR 2 (X — > Y7) ranging from 10 to 60 fb for the general 
case of massive X and Y. For the special case of My pa 0, the 95% CL upper limits on 
a(e + e" -> XX) • BR 2 (X -> Y7) range from 20 to 45 fb, for M x > 45 GeV. These results 
are used to place mo del- dependent lower limits on the Xi mass in a specific light gravitino 
LSP model [5]. Masses between 45 and 99 GeV are excluded at 95% CL. All limits assume 
that particle X decays promptly. 

2 In the 70 — 80 GeV region the limits are actually slightly worse than those along the My=0 axis of 
Figure |21 despite the more efficient background suppression of the M™ ax cut, relative to the kinematic 
consistency cuts applied in the general case. This is due to a deficit of selected events in this region, 
compared to the expected background when using the general kinematic consistency requirements. 
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Sample 


C (pb- 1 ) 


(GeV) 


< Vs > 


N obs 


,^77(7) 

1 *exp 




« 7 2 ? j (pb) 


192 


28.9 


190 - 194 


191.6 


4 


4.26 ±0.11 


0.21 ±0.10 


0.222 ±0.003 


196 


72.3 


194 - 198 


195.6 


5 


9.97 ±0.25 


0.11 ±0.05 


0.215 ±0.002 


200 


74.8 


198 - 201 


199.5 


14 


10.10 ±0.25 


0.29 ±0.08 


0.207 ±0.001 


202 


39.2 


201 - 203 


201.7 


6 


5.21 ±0.14 


0.23 ±0.10 


0.203 ±0.002 


205 


79.1 


203 - 206 


205.0 


10 


10.34 ±0.26 


0.19 ±0.06 


0.198 ±0.001 


207 


132.2 


206 - 209 


206.6 


15 


17.28 ±0.43 


0.17 ±0.04 


0.196 ±0.001 



Table 1: Results of the selection applied to the OPAL 1999 and 2000 data samples. Shown 
for each sub-sample are the integrated luminosity C, the centre- of -mass energy range, the 
luminosity-weighted mean centre- of-mass energy, the numbers of events observed and ex- 
pected, and the measured and predicted cross-section for the process e + e~ — > 1^77(7), 
within the kinematic acceptance of the selection. Predicted values were obtained using the 
KK2f Monte Carlo generator. The errors shown are the sum of the statistical and system- 
atic uncertainties. 



M x (GeV) 


M Y =0 


M Y = M x /2 


M Y = M x - 10 


M Y = M x - 5 


102.5 


74.5 ± 1.2 


74.7 ± 1.1 


63.2 ± 1.3 


33.8 ± 1.5 


100 


74.5 ± 1.2 


74.4 ± 1.1 


61.4 ± 1.3 


32.3 ± 1.5 


90 


74.3 ± 1.2 


75.1 ± 1.1 


60.4 ± 1.4 


36.2 ± 1.5 


80 


73.2 ± 1.2 


73.5 ± 1.2 


65.4 ± 1.3 


37.8 ± 1.5 


70 


74.1 ± 1.2 


71.7 ± 1.2 


62.0 ± 1.4 


39.0 ± 1.5 


60 


73.8 ± 1.1 


71.5 ± 1.2 


62.5 ± 1.4 


41.2 ± 1.5 


50 


72.1 ± 1.2 


71.5 ± 1.2 


65.2 ± 1.3 


43.5 ± 1.5 



Table 2: Selection efficiencies (%) for the process e + e — > XX, X — > Y7 at yfs = 206 GeV 
for various Mx and My (GeV), after application of kinematic- consistency cuts. Not shown 
are the values for M Y = 20 GeV, M Y = M x - 15 GeV and My M x 2.5 GeV. The 
errors shown are due to Monte Carlo statistics only. 





a™"(Mx,M Y ) 


<-(M x ,M Y ) 


< n (Mx) 


a 9 "T(M x ) 


192 


138 fb 


296 fb 


143 fb 


288 fb 


196 


60 fb 


125 fb 


71 fb 


87 fb 


200 


57 fb 


278 fb 


57 fb 


237 fb 


202 


105 fb 


323 fb 


106 fb 


206 fb 


205 


52 fb 


183 fb 


70 fb 


130 fb 


207 


31 fb 


90 fb 


45 fb 


70 fb 



Table 3: Results of individual limit calculations at each centre- of-mass energy. The first 
column shows the data sample. The second and third columns show the maximum and 
minimum 95% CL limits on <r(e + e~ — > XX) • BR 2 (X — > Y7) in the (Mx, My) plane, for 
the case of massive Y for Mx > Mz/2 and Mx — My < 5 GeV. The last two columns 
show the minimum and maximum 95% CL limits obtained for the special case of My ^0, 
for Mx values between 45 GeV and the kinematic limit. 
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M x 


Selection 


Selection efficiency (%) with 


Ndata 




(GeV) 


efficiency (%) 


Mf ax > M x - 5 GeV 






102.5 


75.6 ± 1.1 


73.6 ± 1.3 


2 


1.28 ±0.08 


100 


75.7 ± 1.1 


72.7 ± 1.3 


2 


2.08 ±0.10 


90 


74.9 ± 1.1 


72.5 ± 1.2 


3 


4.14 ±0.16 


80 


73.7 ± 1.2 


71.3 ± 1.2 


4 


6.13 ±0.22 


70 


74.5 ± 1.2 


71.7 ± 1.2 


5 


8.51 ±0.28 


60 


73.9 ± 1.2 


72.2 ± 1.2 


5 


11.25 ±0.34 


50 


72.3 ± 1.2 


69.5 ± 1.2 


10 


14.85 ±0.42 



Table J,.: Selection efficiencies as a function of Mx for the process e + e~ — > XX, X — > Y7, 
/or My ~ ai y^s = 206 GeK T/ie second column shows the efficiency of the general 
selection. The third column shows the efficiency including the additional cut on M^ ax . The 
errors on the efficiencies are statistical only. The fourth column shows the number of events 
from the 205 — 207 GeV data sample consistent with the mass value Mx. The last column 
shows the corresponding number of expected events from the process e + e~ — > ^77(7), 
obtained using KK2f along with the corresponding uncertainty (statistical plus systematic). 
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OPAL 




Figure 1: Kinematic quantities of selected multi-photon events. Shown are a) the recoil- 
mass distribution b) the distribution of the scaled energy of the second photon c) the distri- 
bution of the invariant mass of the 77 system and d) the scaled transverse momentum distri- 
bution for the 77 system. The data points with error bars represent the selected OPAL data 
events. In each case the histogram shows the expected contribution from e + e~ — > ^7777(7) 
events, from KK2f, normalized to the integrated luminosity of the data. The expected back- 
ground from other sources (1.2 ± 0.3 events) is not shown. 
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M x (GeV) 



Figure 2: The shaded areas show 95% CL upper limits on the quantity <r(e + e — > XX) ■ 
BR 2 (X -> Y7) at y/s = 207 GeV obtained from all OPAL data with y/s > 189 GeV, under 
the assumption that the cross-section scales as (3x/s. No limit is set for mass-difference 
values Mx — My < 5 GeV, defined by the lower line above the shaded regions. The upper 
line is for Mx = My ■ 
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Figure 3: The calculated value of M™ x for events selected from a) the 192 — 202 GeV 
data sample and b) the 205 — 207 GeV sample. In each case the data points show the 
OPAL data and the unshaded histogram shows the expected distribution from the Standard 
Model process e + e~ — > 1/7/77(7), evaluated using KK2f and normalized to the integrated 
luminosity of the data sample. In b) the shaded histogram shows the expected distribution 
for the signal process e + e~ — > XX, X — > Y7 for Mx = 100 GeV with arbitrary production 
cross-section. 
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M x (GeV) 

Figure 4: 95% CL upper limits on a(e+e" -> XX) • BR 2 (X -> Y7) at 207 GeF/or My ~ 
obtained from all OPAL data with \/s > 189 GeV. The lightly shaded region shows the 
excluded region obtained using only the OPAL 207 GeV data sample. The darker region 
shows the exclusion region obtained using all OPAL data with y/s > 189 GeV, assuming that 
the cross-section scales as flx/s. The line shows the prediction of an example light gravitino 
LSP model [5]. Within that model, xl masses between 45 and 99 GeV are excluded at 95% 
CL. These limits assume that particle X decays promptly. 
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